Zero-temperature phase diagram of Bose-Fermi gaseous mixtures in optical lattices 
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We study the ground state phase diagram of a mixture of bosonic and fermionic cold atoms 
confined on two- and three-dimensional optical lattices. The coupling between bosonic fluctuations 
and fermionic atoms can be attractive or repulsive and has similarities with electron-phonon coupling 
in crystals. We investigate behavior of the mixtures in the limit, where the Bogoliubov sound velocity 
that dictates bosonic dynamics is comparable to the Fermi velocity, hence the retardation effects are 
important part of the physics. The dynamic Lindhard response function of the fermionic density 
to changes in the bosonic number of particles above some critical frequency can alter the sign and 
in consequence the inter-species interaction between particles becomes repulsive in contrast to the 
static limit (instantaneous and always attractive). Considering the above we show that the structure 
of the phase diagrams crucially depends on the difference in masses of the bosons and fermions. We 
discuss the situations where integrating out fermionic field provides an additional interaction that 
can decrease or increase bosonic coherence. 
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I. INTRODUCTION 



Trapping and cooling Bose-Fermi mixtures of dilute 
quantum gases has opened a wide area of research in 
atomic physics. The interactions between bosonic and 
fermionic species interconnect two systems of fundamen- 
tally different quantum statistics. The diluteness of the 
gaseous mixtures allows one to treat the interactions be- 
tween particles in terms of binary collisions. In conse- 
quence we can replace the real inter-atomic potential by 
a pseudo-potential characterized by only one parameter, 
the s-scattering length. The latter is experimentally[lf 
4 tunable by exploiting optically or magnetically in- 
duced Feshbach resonances[5j. Despite its simplicity the 
interaction potential (mathematically ill defmed|6 ) of 
the ultra-cold multi-component gases confined in opti- 
cal lattices is responsible for a wealth of novel quantum 
phases [7, 8J including charge density waves (CDW)|Slll(J|. 
as well as supersolid behavior |llfH5] . The nature of 
the phase transition and qualitative phase diagram for 
one-component bosonic system can be inferred based on 
very simple arguments [14]. When tunneling between lat- 
tice sites of the bosons is suppressed, compared with 
point-like interaction between them, the system can un- 
dergo a quantum phase transition between a superfluid 
(SF) phase (characterized by large number fluctuations 
at each lattice site), and a Mott insulating (MI) phase 
where each lattice site is occupied by precisely an inte- 
ger number of bosons without any number fluctuations. 
Adding to a bosonic system a fermionic ingredient and 
allowing for the mutual repulsion or attraction between 
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species of different statistics strongly affects the equilib- 
rium properties. Increasing the boson-fermion repulsion 
drives the system towards spatial separation whereas at- 
traction gives rise to implosion[15 16 . The dynamics un- 
derlying the phase transitions in the Bose-Fermi mixtures 
is produced by the small changes of the bosonic density 
which induce a modulation of the fermionic density. As 
a consequence of the feedback of the fermionic pertur- 
bation a shift of the bosonic energy occurs, thereby in- 
ducing an additional attraction or repulsion that changes 
the original interaction between bosons[17] . In the case of 
deep optical lattices and small densities, the coherent de- 
scription of the system provided by the Gross-Pitaevskii 
equation is not reliable due to rising effects of correlation. 
The experimental data clearly demonstrated that adding 
a fermionic cloud to strongly interacting bosons always 
results in a decay of visibility of the interference pattern 
in time-of-flight images [TJ EH H]. Moreover, the scale of 
disappearance of the coherence in the mixtures is dif- 
ferent for attractive and repulsive scattering lengths[4 . 
To predict such behavior theoretically one can include 
the more general, than one-component Bose-Hubbard 
(BH), multi-band model. In mentioned approach if the 
higher-band renormalization of the boson parameters is 
dominant over the fermion screening of the interaction, 
the Mott-insulating lobes in the Bose-Hubbard phase di- 
agram are enhanced for either sign of the Bose-Fermi 
interactions [THJ HSJ. On the other hand inclusions of 
the retardation effects^!]] (which arises from the pres- 
ence of very low energy excitations in a Fermi sea) give 
rise to so-called orthogonality catastrophe[21 . Another 
approach to quantum mixtures of particles of unequal 
masses, when the difference in the tunneling amplitudes 
between heavy bosons and light fermions is large enough 
to neglect quantum nature of the bosons, provides to 



description of the system in Fermi-Bose version of the 
Falicov-Kimball modcl|22_. The mutual interactions of 
bosons and fermions can affect the spectrum of collec- 
tive excitations in the collisionless regime as the mix- 
ture goes toward either demixing or collapse [T5]. It has 
been shown that mode-mode coupling effects may arise 
when sound velocity of the Bose gas is comparable to the 
Fermi velocity of the fermions [23] . The energy spectra 
of the bosonic and fermionic mixtures and phase dia- 
gram were also obtained by the field theory methods |24|. 
however calculated phase boundary does not change the 
structure and only shifts the chemical potential. The 
properly constructed effective theory lead to an effective, 
fermion mediated, long-range interaction between bosons 
with alternating sign that is the origin of the CDW and 
can explain the MI-CD W phase separation [9 . There is 
also another possibility of the analysis of the mixtures of 
atoms with different statistics where the second species 
is strongly localized on random sites which can lead to 
random shifts of the on-site energies and, in consequence, 
the disorder with discrete probability distribution is cre- 
ated [55]. Recently the experiments [26J on a harmon- 
ically trapped mixtures of atomic bose-bose gases show 
that the presence of relevant fraction of the 41 K bosonic 
species modifies the quantum phase transition occurring 
in 87 Rb inducing a significant loss of coherence similarly 
to bose-fermi systems that can be explained in the mean- 
field theory frame work [27]. 

The aim of this work is to study the superfluid to Mott- 
insulator zero-temperature phase transition by means of 
the Bose-Fermi-Hubbard model in two- (2D) and three- 
dimensional (3D) optical condensates. In order to find 
a phase boundary for BF mixtures very sophisticated 
methods and calculations are required. Only several the- 
oretical works concentrated on the inherent difficulty of 
dealing with BF Hubbard Hamiltonian originates from 
the non-perturbative nature of the model and retardation 
effects. To elucidate the quantum phase transition in op- 
tical lattices, where the kinetic energy scale is less than 
the dominating repulsive energy and density-density cou- 
pling between species with different statistics comes into 
play, we have adopted a theoretical approach for strongly 
interacting fermions [SSJ to the BF Hubbard model in a 
way to include the effects of particle number fluctuations 
and make the qualitative phase diagrams more quantita- 
tive [5S]. To facilitate this task, we employ a functional 
integral formulation of the theory that enables to per- 
form functional integration over fields defined on differ- 
ent topologically equivalent classes of the U (1) group, 
i.e., with different winding numbers. An inclusion of the 
winding numbers (comes from periodicity of the phase 
variables) is unavoidable in order to properly construct 
the phase diagram and the Poisson re-summation for- 
mula turns out to be very useful for derivation of the 
topological term of the partition function. The quantum 
rotor representation method we use is deeply rooted in 
the gauge symmetries of the model. We construct an 
invariant theory introducing an appropriate U (1) gauge 



transformation. In Sec. II we review the Hamiltonian 
for the system and show the connections of the parame- 
ters to the experimentally measured quantities in optical 
lattices. Sec. Ill contains description of the method we 
use and can serve a guidance to obtain the critical line 
equation presented in Sec. IV. Before showing the phase 
diagrams for the Bose- Fermi Hubbard model in the quan- 
tum rotor description we make some general remarks in 
Sec. V concern the phase boundary equation and com- 
pare our results with the diagrammatic perturbation ap- 
proach to the one-component Bose-Hubbard model for 
experimentally accessible densities of the particles. Sec. 
VI presents the discussion of the ground state phase di- 
agrams for the Bose-Fermi Hubbard Hamiltonian calcu- 
lated within quantum rotor approach. The Sec. VII. is 
devoted to some concluding remarks. The appendixes 
contain the derivation of relevant formulas of the main 
text and are introduced to keep the text self-contained. 



II. HAMILTONIAN 

For bosons confined in optical lattices the two main en- 
ergy scales are set by the hopping amplitude proportional 
to % (which sets the kinetic energy scale for bosons) due 
to the particles tunneling, and the on-site interaction 
Ub > 0. For if, > Ub the phases of the superfluid order 
parameter on individual lattice sites are well defined. On 
the other hand, for sufficiently large repulsive energy Ub, 
the quantum phase fluctuations lead to complete suppres- 
sion of the long-range phase coherence even at zero tem- 
perature. The competition between the kinetic energy, 
which is gained by delocalizing bosons over lattice sites 
and the repulsive interaction energy, which disfavors hav- 
ing more than one particle at any given site, can be mod- 
eled by quantum Bose-Hubbard Hamiltonian[14 . The 
physics of the bosonic and non-interacting spin-polarized 
(collisions in the s-wave channel are forbidden by their 
statistics) fermionic mixtures with density-density inter- 
action between species Ubf leads to Bose-Fermi-Hubbard 
Hamiltonian [5J: 

- ^2*fijfifj - Vf^n fi + U bf Y^n bi n fi , (1) 

where &](//) and bi(fi) stand for the bosonic 
(fermionic) creation and annihilation operators n^ = 
b\bi, [nfi = fifi) is the boson (fermion) number oper- 
ator on the site i, and the reduced chemical potential 
P-b = Mfe + Ub/2 controls the number of bosons and /// 
fermions respectively. Here, (i,j) identifies summation 
over the nearest-neighbor sites. Furthermore, tuj (tfij) 
is the hopping matrix element for bosons (fermions) . For 
simplicity, we neglect the inhomogeneous magnetic trap 



potential. If the on-site boson-fermion coupling strength 
Ubf becomes very strong the dilute gaseous mixtures are 
unstable to phase separation (Ubf > 0) or to collapse of 
the phase separated configuration (Ubf < 0) |31[3D]. The 
presence of the lattice will introduce kinetic energy scales 
tb(f) competing with Ubf stabilizing the system. We 
assume that an optical lattice created by the counter- 
propagating laser beams is deep enough and we can re- 
strict ourselves to the lowest Bloch bands. The corre- 
sponding experimental parameters can be estimated by 
following relations [12] 
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scattering lengths can be continuously tunc in the 
experiments [THI] inducing attractive or repulsive inter- 
action between species. The k = 27r/A is the wavelength 
of the laser and Ef = H 2 k 2 /2m x is the recoil energy and 
m T is the atomic mass. 



In the next section we will integrate over the fermionic 
fields since the action is quadratic in /j (r) variables. We 
attempt to reduce the large number of degrees of freedom 
in the partition function to the few which dominate the 
low energy physics. 



A. Integration over fermionic fields 

Before integrating out of the fermionic degrees of free- 
dom we write the action in the form: 



S b [b, b, rib] = / dr I ^ &t ( r ) jrb- 
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We notice that adding the inter-species interaction term 
to the fermionic part of the action 



III. DESCRIPTION OF THE METHOD 

We write the partition function of the system switching 
from the particle-number representation to the conjugate 
phase representation of the bosonic degrees of freedom 
using the bosonic and fermionic path-integral over the 
complex fields a, (r) and f+ (r) depending on the "imag- 
inary time" < t < f3 = 1/k^T with T being the tem- 
perature: 



[VbVbVfVf] e 



-Sb.bJJ 



The action S is given by 



where 



S = S B [b,bJ,f]+ / dTH(T), 

Jo 
S B [b,b,fj] = 53/ drh^^-b^r) 



(5) 
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5/+^ = 53 / drdr'fi (t) 



G f+int (r,r') f 3 (r') 



(9) 

allows one to integrate over fermionic fields, because the 
action is Gaussian in fi (r) operators. Resulting partition 
function takes the form 



[VbVb] 



-S b \b,b,n b I Tr In G, 



(10) 



We will be looking for solutions obeying translational 
invariance in the "imaginary time" direction, i.e., such 
that the partition function depends only on the differ- 
ence \t — r'\. Expanding the trace of the logarithm in 
Eq. (flOl we have 



TrlnGT 1 ^ = -Tr InGy - TrG int G f 



1 / - = \ 2 
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G in t(r,T') = U bf bi(r)b s ('T J )6 ij 8{T-'T J ). (12) 

Uj 



We defined I; L j = 1 if i,j are the nearest neighbors and 
equals zero otherwise. Trace over first term of the ex- 
pansion gives constant contribution of the fermions, in 
the non-interacting system, to the action. Second one 
induces a shift in the chemical potential of bosons. The 
third term after exploiting Fourier-Matsubara transform 
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where u t = 2n£/f3 (v t = n{2£+l)/(3) with (I = 
0, ±1, ±2, ...) are Bose(Fermi)-Matsubara frequencies re- 
specting periodic (anti-periodic) boundary conditions of 
the bosonic (fermionic) field operator, reduces to 



TrlnG 
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x x(Q,^.i«/,/3)A_ q (-^), (14) 



where A q (ui) = b q (u>e)b q (u} e ) and \ (q, iv e , p,f,r3) is 
called Lindhard function 



x{q,m,(i>f, 
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(15) 
In the above / [t fk , p f , /3] = 1/ {exp [/3 (t fk ~ (if)] + 1} is 
the Fermi distribution function and tfk is the fermionic 
dispersion relation. To stay in the local (momentum in- 
tegrated) regime we perform q and k integration over the 
first Brillouin zone and in the T40 limit using an ana- 
lytic continuation ivi — >u; + if. we obtain imaginary part 
X" (^,M/) — I m X ( w jA*/) °f the local dynamic Lindhard 
function (see Appendix for details) 

/+oo 
dx [0 (x - J - fi f ) -Q(x- p f )\ 
-co 

x p(x)p(x — u/) (16) 

that satisfies sum rule which is just the conservation 
of the number of particles. In the above p (£) = 
iV -1 J^k $ [£ — tfk] is the density of states and (x) is 
the unit step function. Therefore, the corresponding real 
part x' (Wj Pf) = Rex (w, Pf) can be easily deduced from 
Kramers-Kronig relation 
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(17) 



Later, without any loss of generality, we drop the chem- 
ical potential dependence writing x' ( w j/^/ = 0) = x' ( w ) 
and x"(w,^(y=0) = x" (cj). Finally the "imaginary- 
time" partition function with integrated out fermionic 
degrees of freedom, in the local approximation (see Ap- 
pendix for details), can be written as: 



Z = [VbVb] e 



-S tt\b,b\ 



(18) 
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Figure 1: (Color online) Real tfx' (u) an d imaginary tfx" ( tJ ) 
part of the local dynamic Lindhard function for square (2D) 
and cubic lattice (3D), in the low-temperature limit, as a 
function of frequency, where the normalized fermionic poten- 
tial Pf/tf = is equal zero. The normalized value of the 
critical frequency u -J shows where induced, frequency- 



"crit 

dependent, effective part of the interaction 
changes character from attractive to repulsive. 
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with the effective action 

S eS [b,b] = j drr^bi 
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(19) 



expressed in terms of bosonic degrees of freedom only. 
We want to emphasize that applying a local, in the 
Matsubara-imaginary time, approach we neglect any dis- 
sipation effects. Whereas locality in real space rules out 
some parts of interesting physics such as the charge den- 
sity wave, namely an insulating phase with modulated 
density or the supersolid phase, presenting the coexis- 
tence of superfluidity and a periodic spatial modulation 
of the density, different from that of the lattice. On 



the other hand, the long-range character of the fermion 
mediated interaction between bosons with the fermion- 
induced mean field potential can lead to spatially homo- 
geneous regions of commensurate CDW [5]. The moti- 
vation of the local approximation was the idea that re- 
sponse of an interacting system can be pictured as the 
response of a non-interacting system to an effective self- 
consistent field, that depends on global properties such 
as the particle densities. A question of both fundamental 
and practical interest is, to what extent can the physics of 
the exact non-local interaction be captured by an approx- 
imate local theory? It seems that local approximations 
often work surprisingly well, yielding energies very accu- 
rately, without suffering from some of the characteristic 
drawbacks of non- locality [3T| ■ 

From Eq. ( 19 I it is concluded that there is a striking 



resemblance to the one-component Bose-Hubbard action 
with the original repulsive interaction replaced now by 



U b -> U cS = U b + U^x' (u, Pf) 



(20) 



which is the induced, frequency-dependent, effective in- 
teraction between bosons. From Eq. p0| we see that 
integrating out fermionic field from BF Hubbard Hamil- 
tonian provides an additional interaction among bosons, 
which is not affected by the attractive or repulsive na- 
ture of the inter-species interaction ±£4/. Before we 
proceed with further calculations let us make some re- 
marks. The substitution we introduced in Eq. (20 1 is de- 



ceptively simple and can lead to the assumption that the 
phase diagram of the BFH model can be easily derived 
from the critical line of the one-component BH Hamil- 
tonian, which already has been obtained in several ap- 
proximations. Unfortunately, as we will see in the next 
section, U^x' ( w > Pf) is n °t the only one ingredient to 
the final equation for the critical line. The additional 
part, which influences the phase boundary line condition, 
comes from the inter-species interaction and number of 
fermions added to the system has an impact on amplitude 
of the order parameter. Moreover, the chemical potential 
for fermions fif/tf is shifted (we postpone calculations of 
it now and show proper formula later) because of the in- 
duced effective interaction between them. 

Next step in the calculations depends on the ratio 
mb/mf of the masses of bosons and fermions which can 
be seen from Eq. Q. Consequence of the latter is the 
fact that the speed of the Bogoliubov sound c b for bosons 
differs from the first sound vj of the ideal Fermi gas. In 
typical experimental realizations 40 K- 87 Rb systems the 
acoustic long-wavelength (boson) c b and fermion Vf ve- 
locities are comparable or Cb/vf < 1 (boson mass is 
larger than fermion). Therefore we do not restrict our 
calculations to the static limit but consider also local 
dynamical response function (thus including the retarda- 
tion effects). If c b is much larger than v/ (in the 40 K- 
23 Na species we have c b /vf ~ 5) the resulting interaction 
between bosons is instantaneous and always attractive 
(with U^fX' (Wj/Zf) < 0) so using the static approxima- 
tion is justified with an error which involves the small 
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Figure 2: (Color online) Imaginary tfx"(u,iif) part of the 
local dynamic Lindhard function for cubic lattice in the space 
of the parameters: normalized frequency ui/tj and fermionic 
chemical potential /Lt//t/ in the zero-temperature limit. 
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Figure 3: (Color online) Real tfx' (w ; A 4 /) P ar t of the local dy- 
namic Lindhard function for cubic lattice in the space of the 
parameters: normalized frequency w/t/ and fermionic chem- 
ical potential Hf/tf in the zero-temperature limit. 



parameter Vf/c b . 

The imaginary part of the Lindhard response func- 
tion rises to a broad peak before falling and the real 
part takes zero when changes in x" ( w ) are the biggest 
(see Fig. [I]). The real part x' ( w ) is negative for the 



< 0.4710 = 



,2D 



,3D 



for square and 



frequencies u/8tf 

w/12tf < 0.3525 = 

We normalized the frequency by the width of the band 



for cubic lattice (see Fig. IT]) 



for non-interacting fermions to show the scale of the en- 
ergy. Nevertheless, for higher values of the normalized 
frequency u > w cr i t the induced part of the interaction 
U£fx' (u>) can be positive and increase the effective inter- 
action between bosons, in consequence provide stronger 
localization these species on lattice sites. In the case of 
large fermion hopping and commensurate filling with the 
lattice the effective long-range density-density interaction 
between bosons has alternating sign and is the origin of 
the charge density-wave phases [9]. The higher values of 
the normalized chemical potential for the fermions Hf/tf 
decreases the values of both (Fig. 2] and Fig. [3]) real 
and imaginary part of the local Lindhard function. That 
leads to the situation where the terms containing explic- 
itly the average density of fermions n-p will acquire more 
significance than terms with exclusively the inter-species 
interaction [/&/. 



B. Static and periodic bosonic fields and gauge 
transformation 



Unfortunately the effective action is not quadratic in 
bosonic fields hi and we have to decouple the effective 
interaction term in Eq. ( 19 1 by a Gaussian integration 



over the auxiliary scalar potential fields 



Vi (r) = V l0 + VI (r) 



with static 



and periodic part 
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(21) 
(22) 

(23) 



where uj£ is the Bose-Matsubara frequency. We observe 
now that effective BF Hubbard Hamiltonian has a local 
U (1) gauge symmetry, when expressed in terms of the 
underlying boson variables. This points out a possibility 
of an emergent dynamical U(l) gauge field as a fluctu- 
ating complex field attached to bosonic variables, which 
is dynamically generated, by interacting bosons. Thus, 
the periodic part V( (t) = V{ (r + (3) couples to the lo- 
cal particle number through the Josephson-like relation 
fa(r) = V- (t), where 



h (r) 



_ dfa (t) 
" dr 



= e -Mr) 



ld_ 

i dr 



(24) 



The quantity <f>(r) is the U(l) phase field and satisfies 
the periodicity condition fa (/3) — fa (0) as a consequence 
of the periodic properties of the V( (r) field in Eq. ( 23 1 . 



Next, we perform the local gauge transformation to the 
new bosonic variables 



b(r) 
bi(r) 
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ai (r) 

at (t) 



(25) 



that removes the imaginary term —i Jl drfa (r) nu (r) 
from all the Fourier modes except at zero frequency. 
From the above we deduce bosons have a composite 
nature made of bosonic part a% (r) and attached "flux" 
exp [ifa (t)]. Due to such U(l) gauge invariance, the 
fluctuations and the phase have the dynamics of U (1) 
gauge field. 



C. Gauge group U (1) governed phase only action 

By integrating out the auxiliary static field V^o we cal- 
culate the partition function with an effective action ex- 
pressed in the form of the propagator G 



[Vfa e 



-£ l / "dr[ 5I i-^(T) + i T ^«r)]+'ErlnG- 1 



(26) 

where Hb/Ub — Hb/Ub + 1/2 is the shifted reduced bosonic 
chemical potential. In the above exp I — TrlnG -1 ' ~ 
det G and the determinant takes the form 



dot G 




x a t (r) 



J0i(i") + , .. p -i<t>j(, T ) 



(27) 



We parametrize the boson fields a, (r) = a® + a i (r) and 
incorporate fully our calculations to the phase fluctua- 
tions governed by the gauge group U (1). Assuming non- 
fluctuating amplitude at low temperatures a, (r) = a , 
we drop the corrections, which was proved to be justi- 
fied in the large Ub/tb limit we are interested in |291I32| . 
The amplitude fluctuations are massive one and do not 
play important role in the low energy scales. It is very 
convenient to define the order parameter 

* B = (h (r)) = (a, (t) exp [ifa (r)]) = a ip B , (28) 

which signals the emergence of the superfluid phase and 
vanishes in the Mott-insulator state. The SF state is 
characterized by spontaneously breaking of the U (1) 
symmetry of Bose-Fermi-Hubbard Hamiltonian. Note, 
that a nonzero value of the amplitude ao in Eq. ( 28 1 
is not sufficient for superfluidity. To achieve this, also 
the phase variables <f> in Eq. ( 28 ) , must become stiff and 



coherent, which implies -0b t^ 0. As we see in the next 
sections the presence of the fermions and density-density 
interactions Ubf between species of different statistics can 
also change the amplitude of the order parameter. Af- 
ter mentioned assumption the inverse of the propagator 
becomes 
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T = G7 ) l [1-TGo 



(29) 



The explicit value of the amplitude a® in Eq. (28) 



can be obtained from minimization of the Hamiltonian 
dV. (do) /da = 0. Therefore, we write 
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f = e*** (T) *w<e-^M. 



(30) 
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Expanding the trace of the logarithm in Eq. (26) and 



making use the above we obtain up to the second order 



in the amplitude of the order parameter Eq. ( 28 I 
Tr In G~ x = -Tr In G - TrTG 



(32) 



Trace over first term of the expansion, as previously, 
not containing any fluctuating field variables, gives an 
inessential constant contribution to the action. Let us 
consider the second order term in more detail 





Tr (f G ) = J2 **J / dTdT ' 



where 

<>bij " 



2. 
UQlbij 



Z)(»,i>*Wj+W> u bf 



(33) 



Ub 



Ub 



n, tuj, (34) 



the hopping matrix elements are re-normalized by the 
amplitude of the order parameter. We see that in com- 
parison to pure bosonic case there is an additional shift 
— U b fn-p/U b that depends on the average of the fermion 
concentration and normalized inter-species interaction. 
The above was also obtained in the effective bosonic 
model and recognized as a mean-field contribution [9]. 
Finally, the partition function Eq. ( 26 ) becomes 



Z = / [V<j>] e" 5pha 



(35) 



with an effective action expressed only in the phase fields 
variable 



^phasc 







He 



2U c s i u e s 



t b J2 e 



<t>i( T )j.. e -4>j(T) 



\52(i,j) tbij + Vb - U b fn F ) t b /U b . 



(36) 



The total 



where t b — 

time derivative Berry phase imaginary term in Eq. ( 36 1 is 



nonzero due to topological phase field configurations with 
4>i(l3) — 4>i (0) = l-Kwii [rrii — 0,±1,±2...) that results 
in topological ingredients to the correlator we will see 
below. Therefore, we concentrate on closed paths in the 



"imaginary time" (0, /3) labeled by the integer winding 
numbers m^. The path-integral 



[V4>] 



E 



2;r 



[2ty(0)] 



(j)i(r)-\-27Tmi 



pty(r)]. 



(37) 
includes a summation over rrii and in each topological sec- 
tor the integration goes over the gauge potentials. There- 
fore, we do not ignore the compactness of the gauge fields. 
To proceed, we replace the phase degrees of freedom by 
the uni-modular scalar complex field -0 which satisfies the 
quantum periodic boundary condition ipi (/?) = ipi (0) . 
This can be conveniently done using the Fadeev-Popov 
method with Dirac delta functional resolution of unity, 
where we take tp as a continuous but constrained (on 
the average) variable to have the uni-modular value. We 
introduce 



i = yv^w(z>( 

x 5 (i>i - e**< r >) S (V* - 




(38) 



and 



*(X>( 7 



-N = 



1 

2~7ri 



-\-ioo 



d\ 



x e 



/<f*rA(X\hMr)| 2 -iV) 



(39) 



where N is the number of lattice sites. Introducing the 
Lagrange multiplier A, which adds the quadratic terms 
(in the ip fields) to the action we can solve for the con- 
straint. The partition function can be rewritten to the 
form 
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2iri 



-\-ioo 



e- XN dX 



[VipVip* 



exp < - 5^ / drdr'-ipi [{ibhj + A%) 8 (i 
( id Jo 



hAt,t')}tP*} 



t') 

(40) 



where 



7tf (r, t') = (exp {-i [fr (r) - ^ (r')]}> (41) 



is the two-point phase correlator associated with the or- 
der parameter field, where (• • • ) denotes averaging with 



respect to the action in Eq. (36 1. Because the values 



of the phases <j> which differ by 2ir are equivalent we de- 
compose phase field in terms of a periodic field and term 
linear in r: 



4>i (r) = tpi (r) -I- -^miT 



(42) 



with 4>i (/3) = fa (0) . As a result the phase correlator 
factorizes as the product of a topological term depending 
on the integers rrii an d non-topological one: 



7«(r,r')=75(r,r')7iT(r,r') 



(43) 



Performing the Poisson re-summation formula in 



^ In- n-'\ M m *I 6 

lii ( T ' T ) = — 



^(T-r'W-^E,[t£ ? ' n ?+f&»w] 



E[mj] 



/ii ^— 'I U„ff % % V„cc 



(44) 
and the functional integration over the phase variables 

lij ( r , T ) = 



/PVle-S'^tf^M 



the final formula of the correlator takes the form 



(45) 



7«(t,t') = 



^( 7r £ f - + 7r ^> e " ctf " 



Ueff ' 



x exp 



t/eff 



/8 



,(46) 



where ■& (z, q) is the Jacobi theta function, which comes 
from the topological contribution - summation over in- 
teger winding numbers. The function $ (z, q) is defined 

by 



+oo 



${z,q) = 1 + 2^2 cos ( 2nz ) < 



(47) 



n=l 



and is /3-periodic in the "imaginary time" as well in the 
variable ftb/U e g with the period of unity which empha- 
sizes the special role of its integer values. After Fourier 
transforming one obtains 



tbk = 2tb (cosk x + cos k y ) and three-dimensional tfj^ — 
2t\, (cos k x + cos k y + cos k z ) lattice. 

IV. CRITICAL LINE 

Within the phase coherent state the order parameter 
i{>B is evaluated in the thermodynamic limit N — > oo by 
the saddle point method 5 J- '/ '<5A = and the uni-modular 
condition of the U (1) phase variables translates into the 
equation 



i-i%= lim _^r fc M, 



with 



r -i 

1 k 



jv->oo 7V/3 






(51) 



(52) 



The phase boundary is determined by the divergence of 
the order parameter susceptibility Tk^o {w£=o) = 



Ao-Cfe=o + 7- 1 (^=o) = 



(53) 



which determines the critical value of the Lagrange pa- 
rameter A = Ao and stays constant in the whole global 
coherent phase. To proceed, it is desirable to introduce 
the density of states 



■®4?'(<-£) 



(54) 



because the analytical expressions we use can be advan- 
tageous in evaluating sums over momenta. The corre- 
sponding formulas for square lattice can be written as 



hj K) 



where 



1 4 



E— , 

M 1 " 4 [Ei 



*^(-+er 



U e ff (Jeff 



Z = £e 



-^ £*("»«+& 



(48) 
(49) 



is the partition function for the set of quantum rotors. 
The action Eq. (36), with the topological contribution 



Eq. ( 48 1 , after Fourier transform, is written as 



Scff [ip, ^]=^X! ^ k ^ r k 1 ( W£ ) ^ ( Wf ) ' ( 50 ) 



where T^ (cJe) = A — i&k + 7 _1 (w^) is the in- 
verse of the propagator and £&& is the Fourier trans- 
form of the bosonic hopping matrix elements for two- 
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(0 = 



1 
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4*7, 



and for simple cubic geometry takes form 



(55) 



„3D 



(0 = 




(56) 



with ai — min (— 1, — 2 — £/2£(,) and a 2 = 
max(l, 2 — £/2tfc); K (x) is the elliptic function of the 
first kind. [33 . After summation over Bose-Matsubara 
frequency and for zero temperature limit (3 — > oo we 
can rewrite the critical line equation to the form that 
represents solution of the BF Hubbard model in terms 
of re-normalized pure Bose-Hubbard Hamiltonian in the 
quantum rotor approach: 



v-i 



+00 



piOdt; 



Jnt™-o(2zb + ft-r,+h)hk+*(hft) 



(57) 



In the above v (nb/atUb) = frac (nb/aUb) — 1/2, where 
frac (x) = x — [x] is the fractional part of the number and 
[x] is the floor function which gives the greatest integer 
less than or equal to x\ £ m ax stands for the maximum 
value of the bosonic dispersion spectrum tbk and z is 
the lattice coordination number. The renormalization 
parameters are defined as: 



0.07 



a 



l + ^x'^) 



Ubf 

u b 



"F 



(58) 
(59) 



and allow us to see how adding free fermions to strongly 
interacting bosons confined in optical lattice influences 
the phase boundary. 



V. PHASE DIAGRAMS - BH MODEL 

The zero-temperature phase diagram of the Bose- 
Fermi-Hubbard model Eq. (TTT) can be calculated from Eq. 
( |57| and usually is plotted as a function oitb/Ub, with the 
density of the bosons controlled by a chemical potential 
Hb/Ub- The presence of the fermions implicates two ad- 
ditional different parameters that can by varied namely 
a and 77 in Eq. ( 58 I and Eq. ( 59 ) . The strength of the 



inter-species interaction influences both of them, however 
the sign of Ubf and the average density of fermions n-p 
added to the system affects only 77. Besides, in the local 
dynamic approach the sign of the density-density interac- 
tion depends also on the normalized frequency oj/tj . For 
a general choice of parameters, Eq. (57 1 is easy to solve, 



however considerations of special cases can provide more 
insights into the solution of the problem. In discussion 
we will follow the scheme 



a < 1 and 



Cb/vf > 1 for uj = 

Cb/Vf ~ 1 for W < Wcrit 



77 ^ (60) 

a > 1 and Cb/vf ~ 1 for u > w C rit 

firstly choosing the sign of the 77 and later a in the static 
(w = 0) or dynamic (u ^ 0) limit (see Fig. HI. 

Before we proceed with analysis let us introduce the 
notation for the maximum of the critical value for pa- 
rameter tb/Ub (as a function of the normalized chemical 
potential Hb/Ub) at the tip of the nth (77,3 — n) MI lobe 
for different lattice geometries and model parameters a 
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Figure 4: (Color online) Phase diagrams (tb/Ub-Pb/Ub) f° r the 
square (2D) lattice for different a — 0.5 (higher panel, u < 
w C rit)i a = 1-5 (w > oj cr it) and r\ = —2 (negative scattering 
length a,bf < 0). Dashed line stands for the phase boundary 
of one-component Bose-Hubbard model. Within the lobes the 
Mott insulator phase takes place with *I/b = 0. 



and 77 as follows 



x n (a, 77) = max 



2D, 3D 



(61) 



cntJ a,r] 



The above determines when the transition from MI to 
SF occurs. Values a = 1 and 77 = stand for the one- 
component bosonic case. In Table [I] we show compari- 
son of x n (1,0) for higher densities of the particles calcu- 
lated in the quantum rotor approach (QRA) to very ac- 
curate, recently developed, diagrammatic perturbation 
approach|34| to Bose-Hubbard Hamiltonian. The re- 
sults for 3D BH model obtained in both theories are 
very close and also comparison to quantum Monte-Carlo 
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Figure 5: (Color online) Phase diagrams (tb/Ub-fJ,b/Ub) for 
the cubic (3D) lattice for different a = 0.5 (higher panel, oj < 
k'crit), a — 1.5 (oj > aj cr it) and n = —2 (negative scattering 
length 06/ < 0). Dashed line stands for the phase boundary 
of one-component Bose-Hubbard model. Within the lobes the 
Mott insulator phase takes place with $b = 0. 



(QMC) numerical calculations[35 indicates that methods 
we use are able to properly catch the interesting physics 
of strongly interacting systems. However, we want to 
analyze the phase boundary for number of particles per 
lattice sites higher than one ne > 1 that adds another 
dimension to the analysis and is difficult for the QMC to 
catch. The phase boundary for square lattice shows that 
QRA works well also in low-dimensional geometries, es- 
pecially for higher densities. Nevertheless, the structure 
of Eq. (57 1 can cause some problems when both a and r\ 



are nonzero and the number of bosons is equal one per 
lattice sites. We expect that for r\ > 1 some artificial 
effects may arise for values of the normalized chemical 
potential (ib/Ub ~ close to zero. 



VI. BOSE-FERMI-HUBBARD PHASE 
DIAGRAM 

In the experiments|4j for a degenerate mixtures of 
4 x 10 5 87 Rb bosons and 3 x 10 5 40 K fermions the 
scattering length a b f (and in consequence interaction 
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QRA 
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0.0136102 


0.0103755 


0.0042086 



Table I: Comparison of the maximum of the critical value for 
(tb/Ub) CTit (as a function of the normalized bosonic chemical 
potential fib/Ub) at the tip of the nth (ns = 1t4 and 10) 
Mott insulator lobe for the square (2D) and cubic (3D) lattice 
in the one-component Bose-Hubbard model: DPT - diagram- 
matic perturbation theory [34 ), QRA - our calculations using 
quantum rotor approach). 
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Figure 6: (Color online) Maximum value of the (ib/^ r b) max 
for different a and n = 0.001 (very small amount of fermions) 
for square (2D) and cubic (3D) lattice with n-g. = 1. Vertical 
solid line stands for the {tb/Ub) mgjc obtained in one-component 
Bose-Hubbard model with one particle per lattice site. Above 
the curves the superfluid phase takes place with $b 7^ 0. 



Ubf, see Eq. (Wl)) can be continuously tune between 
— 170ao -j- +800ao below and between — 800ao -= — 200ao 
above Feshbach resonance, where a,Q is the Bohr radius. 
The form of the parameters we choose Eq. ( 58 I and Eq. 



(59 1 allows for its interpretation. The periodicity of the 



phase diagram can be easily deduced from the periodic 
properties of the propagator Eq. (48 1 and strongly de- 



pends on a however, the interaction between species does 
not generate additional Mott lobes in the phase diagram. 
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The above is in contrast to the strong coupling expan- 
sion and exact diagonalization method applied to the 
system of two alkali-metal atoms with different masses 
where (for very small lattice sizes and quenched disor- 
der) the MI phases with integer filing factors disappear 
for boson-impurity interaction energy larger than on-site 
atom-atom interaction energy itself and also the MI phase 
exists for incommensurate bosonic filling [25 . If we fix 
the number of fermions n-p and inter-species interaction 
Ubf in Eq. ( 59 ) still there is a dynamic part of the lo- 



cal Lindhard function x' ( W )M/) we have to take into ac- 
count. In the static limit uj/tf — > (where the Lindhard 
response function is purely real) there is nothing unex- 
pected in the behavior of the critical line (see discussion 
below). However, we must stress that even we left the 
frequency dependence apart, there is still very interest- 
ing part of physics remained, because the Lindhard re- 
sponse function for the system with regular density of 
states shows logarithmic divergence as temperature goes 
to zero. These singularities give rise to instabilities in 
the system towards two new ground states a phase sep- 
arated state or a supersolid phase |11II36| . On the other 
hand the oscillation of the induced effective interaction 
between bosons is the origin of the formation of charge 
density waves [9]. 

Taking cj < w cr i t we recover the previous theoretical 
results where, after adding fermions to the system, the 
effective interaction U e g becomes smaller than repulsive 
energy Ub for bosons only (see Fig. El and Fig. |5| and 
superfluid phase increases. The above is best shown on 
Fig. [6] where for very small amount of fermions rip the 
parameter a < 1 causes decreasing the Mott insulator 
region of the phase diagram in comparison to the pure 
bosonic case. However, in the local dynamic limit, when 
u > w cr it the Mott insulator phase is becoming stronger 
and bosons tend to localize on the lattice sites in both 
2D and 3D cases (see Fig. [6] Fig. [7] and Fig. |j. 

One may argue that the Lindhard response function 
depends not only on the frequency u/tt but also on the 
chemical potential fXf/tf for fermions and so far we did 
not restrict ourselves to any particular value of it. In 
many approaches it is a little tricky to handle because 
despite of the absence of any direct interaction between 
fermions [/// = the density-density fluctuations can in- 
deed induce some effective interaction between fermionic 
species [37]. Therefore problem becomes complex and 
many theories just take half-filled band with 7i F = 1 so 
that fif /tf = 0. Alternative approach comes from partial 
particle-hole symmetry Hamiltonian Eq. (flj possesses. 
To make our approach self-consistent we can calculate 
how does a particular value of the fermionic chemical po- 
tential change in the effectively interacting system. We 
remind that the amplitude of the order parameter was 
obtained from minimization condition, assuming nonfluc- 
tuating bosonic amplitude at low temperatures. By op- 
erating a similar procedure we get a shift of the chemical 




0.14 , 



2.5 /v 3.0 

r b b 




0.0 0.5 1.0 1.5 2.0 2.5 u/T j3.0 



h /U b ' 



Figure 7: (Color online) Phase diagrams (tb/Ub-^b/Ub) for 
the square (2D) lattice for different a — 0.5 (higher panel, 
lu < tJcrit), a = 1.5 (ui > aj cr it) and 77 = 1 (positive scattering 
length a,bf > 0). Dashed line stands for the phase boundary 
of one-component Bose-Hubbard model. Within the lobes the 
Mott insulator phase takes place with <1/b = 0. 



potential for fermions 



fXf-^Hf 



U bf I 2z 



tb ^b 

'Ub Ub 



n- 



(62) 



that in the non-interacting case Ubf — reduces to that 
of free particles obeying fermionic commutation relations 
(see also Appendix). 

There is the limit where the system containing gaseous 
mixtures has the same value of x(a, 77) as in the case 
of only bosons confined in optical lattice (see Fig. I9J. 
Again, we take advantage of the choice of the parame- 
ters, that suits well our goal, and make notation of the 
condition very simple. If r\ = 1 — a we have 



X n (1,0) = x n (a, I - a) . 



(63) 



in terms of the original variables, leads to 
2 ^' ''■' " -^ /Ub- For cubic lattice the for- 



The above 

U bf n F /U b = UgfX 1 (u, Vf) I u b 

mula Eq. (63) seems to not hold (Fig. [9l however is 
accurate with a numerical error less than 0.17 percent. 
Therefore, if the number of fermions added to the system 
is equal to the inter-species interaction then bosons be- 
haves as if were unaffected by the presence of fermions. 
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Figure 8: (Color online) Phase diagrams (tb/Ub-^b/Ub) for 
the cubic (3D) lattice for different a = 0.5 (higher panel, 
uj < aj cr it), a — 1.5 (w > o; cr it) and 7/ = 1 (positive scattering 
length a,;,/ > 0). Dashed line stands for the phase boundary 
of one-component Bose-Hubbard model. Within the lobes the 
Mott insulator phase takes place with 4<b = 0. 
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Figure 9: (Color online) The maximum of the critical value 
for the parameter {tb/Ub) mlal with jib = 1 for different a — 4, 
3, 2, 1, 0.75, 0.5, 0.25, 0.1 for square (2D) and cubic (3D) 
lattice. Vertical solid line stands for the (tb/Ub)„ obtained 
in the one-component Bose-Hubbard model with one particle 
per lattice site. Above the curves the superfluid phase takes 
place with ^b t^ 0. 



As a matter of fact we have to remember about sign of 
the scattering length a?,/ and normalized frequency w/£/ 
that also modifies the introduced condition. The recent 
experiment [4] shows that there is an asymmetry in pro- 
files of visibility of the interference pattern (recorded by 
absorption imaging) versus the inter-species scattering 
length that increases with lattice depth. Presented data 
indicate that visibility shows a maximum at the position 
consistent with a b f — 0. Besides, there is an asymmetry 
in a shift of the MI to SF transition boundary. Our cal- 
culations can reproduce latter however if the sign of the 
inter-species interaction is negative the MI phase dimin- 
ishes and quite oppositely for positive scattering length 
and, as we expected, some anomaly appears at the point 
with Hb/Ub = (see Fig. [7]and Fig. [8]). There is no phys- 
ical reason for the phase boundary to change a position 
where the chemical potential for bosons is zero and the 
repulsive interactions are very strong Uj, — > oo. More- 
over, in that case the value obtained from Eq. (57 1 at 
mentioned point (/ib/C/j, = and fixed r\ > 0) is con- 
stant in whole nonzero range of the parameter a and de- 
pends only on the considered topology of the system. The 



similar to fermion-boson loss of coherence for the boson- 
boson species was found using the Gutzwillcr mean-field 
approach |27 . The main effect of the addition different 
species of the same statistics is that the new structure of 
wedding cake appear but the oscillatory behavior of the 
relevant condensate fraction does not necessarily result 
in increase of the overall coherence of other species. The 
later is limited exclusively to the shallow lattice depth 
and was never observed in the experiments. 

We want to stress that one have to be careful with 
the analysis of the phase diagrams. The summary of our 
results for square lattice (we omit qualitatively similar 
results for cubic geometry) is presented on Fig (10 1 and 
Fig (111. For the static and dynamic limit, but below 



the critical frequency to < uj cl -it the Mott insulator region 
on the phase diagram broadens only when the scattering 
length is positive (the part of the surface above the plane 
of the critical value of X\ (1, 0) for the one-component BH 
model). When we must take into account the difference 
in the inter-species masses irib/mf ^ 1 the sign of the 
real part of the local Lindhard response function may 
become positive and, in consequence, a parameter takes 
values above one Fig. (Ill leading to higher repulsive 
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Figure 10: (Color online) The maximum of the critical value 
for the parameter (W-M ma x with the negative real part of 
the local dynamic Lindhard function \' ( aj ) < ( w < Went), 
for cubic (3D) lattice in the space of the dimensionless pa- 
rameters a — rj with ub — 1. The flat surface stands for 



x\ (1,0) (see also Eq. 63 1 in case of the one-component Bose- 



Hubbard model with one particle per lattice site. The dashed 
line stands for the condition where the system of gaseous mix- 
tures has the same value (tb/Ub) mgx as only bosons confined 
in optical lattice. Above the surfaces the superfiuid phase 
takes place with *&b ^= 0. 



Figure 11: (Color online) The maximum of the critical value 
for the parameter (tb/Ub) with the positive value of the 
real part of the local dynamic Lindhard function \ ( w ) > 
(uj > Worit), for cubic (3D) lattice in the space of the di- 
mensionless parameters a — rj with ub — 1- The fiat surface 
stands for a;? 13 (1,0) (see also Eq. 631 in case of the one- 



component Bose-Hubbard model with one particle per lattice 
site. The dashed line stands for the condition where the sys- 
tem of gaseous mixtures has the same value {pb/Ub) as 
only bosons confined in optical lattice. Above the surfaces 
the superfiuid phase takes place with ^b t^ 0. 



energy between bosons even if measured scattering length 
is negative. In that case for a > 1 — rj there is always a 
shift for higher values for the parameter x n {a, rj > 1 — a) 
results in stronger localization of the bosons after adding 
fermions to the system. 



VII. CONCLUSIONS 

It is well known that the ground state of a system of 
repulsively interacting bosons confined in a periodic po- 
tential can be either in a superfiuid or in a Mott-insulting 
state, characterized by integer boson densities. Because 
the phase of the order parameter and the particle num- 
ber, as conjugate variables, are subject to the uncertainty 
principle A0An ~ h, so the bosons can either be in the 
eigenstate of particle number or phase. The eigenstate 
of phase is a superfiuid and that of particle number is 
a localized Mott insulator. Therefore, the quantum MI- 
SF phase transition takes place as the particle density is 
shifted thus facilitating emergence of the superfiuid from 
the Mott insulating state. Adding to bosons particles of 
different statistics and allowing for the mutual repulsion 



or attraction between species strongly affects the equi- 
librium properties. We presented a field-theoretic study 
of the ground-phase diagram in quantum two- and three- 
dimensional gaseous Bose-Fermi condensates where men- 
tioned emulation takes place. We calculated the phase 
diagram using the quantum rotor approach that can re- 
produce the asymmetry in a shift of the MI to SF tran- 
sition boundary for positive and negative inter-species 
scattering length. Analysis of the local dynamic Lind- 
hard function revealed the critical value of the frequency 
for the collective excitations, where the real part of the 
response function (and in consequence the interaction 
between bosons and fermions) alters sign. The choice 
of the parameters of the model led to simple condition 
for the experimentally accessible parameters within the 
phase diagram for Bose-Fermi mixtures is qualitatively 
the same as for one-component repulsively interacting 
Bose system. We also compared the maximum of the 
critical value for tb/Uj, parameter (as a function of the 
normalized chemical potential jib/Ub) at the tip of the 
nth MI lobe for square and cubic lattice with numerical 
diagrammatic method and found them in a good agree- 
ment especially for higher, experimentally realizable, fill- 
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ing factors. The nice feature of presented approach, de- 
scribed in details above, is that all the expressions and 
handling are analytic. It is also worth to notice that 
provided local approximation can be very useful in vari- 
ous situations whenever the retardation effects has to be 
taken into account and we are not interested in effects 
caused by non-locality. 



Appendix A: Local dynamical approach 



The third term of the trace (Eq. 14 1 can be written 
after Fourier transform in form 



Tr\n G]l int = £ 



J k-k 



(ui - ue>)b k _ k > (ut - ui>) 



hh'lj 



tfk + ivt 



frfc-fc' (Wg - m) b k-k' {ui> - ui) 
tfk' + ivt> 

i 

^2 A q (u e ) A_ q (-u e ) . 



(Al) 



In the above we picked up some special value of the fre- 
quency V(. Now, doing the inverse Fourier transform and 
using gradient expansion 

bj (r) = h (r) + (r - r') d T b (r) + O f(r - t') 2 ] (A2) 



we obtain local, in the Matsubara-imaginary time, 
quadratic form of the trace in the bosonic variables 



Tr In G f l int -^(ii/^p/, 



■£ 



dT[6i(r)6i(r)] 2 . 

(A3) 
We performed an expansion were not the bj (r') degree of 
freedom itself but rather its gradients of d T bj (r) are as- 
sumed to be small. The explicit formula of the imaginary 
part for the dynamical Lindhard function is calculated in 
the next section. 



Appendix B: Local (momentum integrated) 
Lindhard function 



To stay in the local regime and using an analytic con- 
tinuation ivg — Hj + ie , where it comes from the causality 
relation of the response function 



lim 



1 



»o+ u ± ie 



= V 



± iirS (w) , 



(Bl) 



where symbol V denotes the Cauchy principal value 
which prevents divergence when u = u/, we calculate 
the explicit value of the imaginary part that is somewhat 
easier to obtain than the real part 



x" (ive^oj + ie,n f ,(3) = ^ 

U.q 



f[tfh,tif,0]~f[tfk+q,Hf,0] 



lim 



tfk — tfk+q — W + ie 



>o+ {2-Kf J_„ 



tfk — tfk+q — W + ie 



.in, r^r^^j^-^^p^piy) 



+0+ 

-oo W - u 

I 



x — y — u + ie 



(B2) 



In the low temperature limit T — » the Fermi distribu- that satisfies sum rule f_ dujA(u) = 1. Therefore, we 



tion becomes / (x) = 1 — (a;) and we write 



can also calculate the real part 



A{u',Hf) = *'»/) 

/+oo 
dx [0 (x — p/ — u>') — (a; — /it/)] 
-oo 

x p (x) p (x - w') (B3) 



x'(w,p/) = 7? 



+°° ^, ,' 



2 

7T./,, 



dw' 4 (u/, fj,f) 

oo 1" W'-W 

+oo ,/-,// I 1 , ,/ 



^^yW . (B4) 



From these results one finds that x'i^^f) 1S an even 
function of frequency while \" { u i P/) i s odd. 
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Appendix C: Fermionic number of particles in the 
effectively interacting system 

In the non-interacting case the number of fcrmions on 
the lattice can be calculated as follows 



with p (£) being the density of states for chosen lattice 
geometry. Introducing a shift Eq. ( 62 1 in the above we 



are able to obtain the particular value of the chemical 
potential for fermions with effective interaction induced 
by the coupling with bosonic species. 



"F 






1 



i\r^exp[0(t /fc -M/)] + l 

„ 1 r+°° 

T->0 1 



/+oo 
«p(0 ©(?-/*/)• (ci) 

-oo 
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